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I review the current status of lattice calculations of the properties of bound states containing one or more 
heavy quarks. Many of my remarks focus on the heavy-light leptonic decay constants, such as /b, for which 
the systematic errors have by now been quite well studied. I also discuss B-parameters, semileptonic form 
factors, and the heavy-light and heavy-heavy spectra. Some of my "world averages" are: fs = 200(30) MeV, 

f B \/%^ = 230(40) MeV, f B Jf B = 1.16(4) and f B , y/%Tjf B \f&7 d = 1.16(5). 



1. INTRODUCTION 

Many of the parameters of the Standard Model 
can be constrained by measurements of the prop- 
erties of hadrons containing heavy quarks. To 
take advantage of such experiments, however, 
one needs theoretical determinations of the cor- 
responding strong-interaction matrix elements. 
Lattice gauge theory provides, at least in prin- 
ciple, a means of computing hadronic matrix ele- 
ments with control over all sources of systematic 
error. Here, I review the current status of these 
computations. 

A review like this is necessarily somewhat id- 
iosyncratic in the topics it covers and the time 
spent on each. I devote a large fraction of my 
time to Jb- Because lattice data for leptonic de- 
cay constants is more extensive than for any other 
heavy quark quantity, I am able to discuss in de- 
tail several vexing issues: renormalon effects, sys- 
tematic errors due to large lattice quark masses, 
extrapolation to the continuum, and quenching 
errors. 1 then turn to the "B-parameters," par- 
ticularly Bs d and Bg a , and to semileptonic form 
factors for B — ► 7r. Finally, I discuss spectra of 
heavy-light and heavy-heavy hadrons. I conclude 
with a few remarks on the unitarity triangle. 

Heavy quark masses, which are included in Lu- 
bicz's talk at this conference pL are examined 
here only for the light they shed on renormalon 
effects. The heavy quark potential is discussed 
in detail by Bali and I omit it here. Re- 
cent reviews of heavy quark physics on the lattice 
by Hashimoto || and Draper Q] complement the 
current treatment. 



2. LEPTONIC DECAY CONSTANTS 

Table [j] and Fig. [l] show recent results for fs ; 
while Fig. || plots the quenched Jb data as a func- 
tion of lattice spacing a. One notices immedi- 
ately that the UKQCD HQ and CP-PACS re- 
sults 0,|l^l (with heavy clover quarks and with 
NRQCD) are somewhat high compared to those 
of the other groups. Furthermore, the errors 
given by CP-PACS are rather small. Does this 
mean one should increase the central value and 
lower the error significantly from the "world av- 
erages" for quenched fg quoted recently by and 
Hashimoto [| (f B = 170(20) MeV) or Draper g 
(f B = 165(20) MeV)? 

In the case of UKQCD, a major cause of the 
quite high result is the choice of the ro parame- 
ter jl5j from the static potential to set the scale. 
Since aro can be extracted easily and precisely 
from lattice data, it provides an excellent way to 
check the scaling of physical quantities as a is 
changed. However, I would argue that because 
ro is not directly determined by experiment, but 
only through phenomenology, aro should not be 
used to set the absolute scale of a. Indeed, Som- 
mer [^5| remarks that an error of roughly 10% 
should be associated with his value ro = 0.5 fm. 

Of course, in the quenched approximation, the 
scale set by different quantities will differ. But 
the experimental uncertainty in ro introduces an 
unnecessary additional error. Indeed, while scales 
set by various directly measurable quantities, i.e., 
f n Q, m p @, /k and m K Q , or m K » and 
m,K jl2|], differ by at most 7% at /? = 6.2, the r 
scale differs from these by as much as 15%. The 
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Table 1 

/b and f B J f B from various groups, 
others use relativistic formalisms. A * 



Computations using NRQCD are explicitly noted with (NR); 
denotes a preliminary result. 



all 



UKQCD result f B = 218(5)+^ MeV || becomes 
195(8)1^ MeV Jl| when f„ is used to set the 
scale. Other systematic errors (some of which 
are discussed below) then easily account for any 
remaining difference between this result and those 
of other groups. 

Turning now to the CP-PACS results 
we see from Fig. ^ that these are from rather 
coarse lattices (a>0.57(GeV)^ 1 ) and that the lat- 
tice spacing dependence is quite large. Although 
the Iwasaki gauge action [^2| used here generally 
gives quite small scaling violations, decay con- 
stants are an exception |23J. In particular, / T has 
an a-dependence that is qualitatively very similar 
to that of f B . Indeed, in the heavy-clover case, 
/b / f-n scales considerably better than /s / m p and 
gives f B = 174(4) (3) MeV instead of the quoted 
188(3)(9) MeV. Similarly, with NRQCD, f B /U 
gives f B = 169(7) MeV instead of the quoted 
191 (5) (11) MeV. In the absence of data at smaller 
a, I believe the /b/Ztt results are more reliable. 
Furthermore, the discrepancy suggests that the 
systematic error estimates should be increased. 



The CP-PACS results are then completely com- 
patible with the average from other groups and do 
not imply a significant increase in central value or 
decrease in error for quenched f B - 

Before determining a world average for f B , I 
now examine several other sources of systematic 
error that affect calculations by various groups. 

2.1. Renormalon Shadows 

In lattice HQET NRQCD ||, and the 

Fermilab |2(| formalisms, corrections to the static 
limit require the addition of higher dimensional 
operators, with coefficients that depend on the 
heavy quark mass. Such power corrections are 
most easily discussed in HQET. The framework I 
will use is that of Martinelli and Sachrajd a ]27| ]; 
for additional relevant discussion see Ref. ]28||. 

In lattice HQET a physical quantity "P, which 
depends on the heavy quark mass mg, can be 
expressed through order 1/mg as 

V(m Q ) = C 1 {\0 1 \) a + ^-(\0 2 \) a (1) 

where the short distance coefficients C\ and C2 
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Figure 1. Results for fs from various groups in 
the quenched approximation (circles) and with 
Nf — 2 (crosses). Solid lines show statistical er- 
rors; dashed lines show systematic and statistical 
errors, added in quadrature. A * denotes a pre- 
liminary result. Calculations using NRQCD are 
marked (NR); all other computations using rela- 
tivistic formalisms. For references see Table @. 



arc pcrturbatively calculable functions of arriQ 
and the coupling g. The higher dimension op- 
erator 02 mixes with 0\ with a power divergent 
(like 1/a) coefficient. In particular, we have 



Ci(am Q ,g) = a(am Q ,g) + 



Ci(am Q ,g) 



ci{am Q ,g) 
ci(airiQ,g) 
C 2 (am Q ,g) 



2amQ 
1 + #.9 2 + + • • • 
#.9 2 + + • • • 
1 + #.9 2 + #.9 4 + • ■ • , 



(2) 



where stands for logs of arriQ or constants. 
For future reference, I denote C2O2/ '(2m<g) as the 
"bare" 1/rag operator. Terms in C\0\ generated 
by loop corrections to the bare 1/m.Q operator I 
call the "subtraction." (The subtraction is just 
[c 1 /(2am Q )]0 1 in the HQET case.) Finally, the 
sum of the bare 1 /rriQ operator and the subtrac- 
tion gives the "renormalized" 1/toq operator. 

It is widely believed, but not rigorously proven, 
that the series c\ and c\ have renormalon ambi- 
guities at C(Aqcd/"iq)- In the sum Ci, these 
ambiguities should cancel. However, since the 
cancellation only occurs at high order, one might 
anticipate p7[ that the low order series for C\ 
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Figure 2. Recent world data for fg in the 
quenched approximation vs. lattice spacing. Only 
statistical errors are shown. For references see Ta- 
ble n 



"converges" rather slowly. I call this effect the 
"renormalon shadow:" although the renormalons 
are formally gone, their influence lingers. We ex- 
pect renormalon shadows in any lattice calcula- 
tion where power law divergences are subtracted 
perturbatively. 

As an example, consider the lattice HQET 
computation of the b quark mass. In the static 
limit, mi, is just the meson mass Mb- A non- 
trivial calculation must therefore include the 
1/mg correction, i.e., the term of O(Aqcd)- The 
relevant coefficients were calculated to two loops 
and combined with lattice data in |2{| • The "bad 
news" from this calculation is that the pertur- 
bative error on nib that would result from stop- 
ping at one loop is ^300 MeV, namely the size of 
the effect one is trying to compute. The rather 
large error can be taken as evidence that the 
renormalon shadow is real in this case. How- 
ever, I also find "good news" in Ref. ^9|: The 
~ 300 MeV one-loop error could have been cor- 
rectly estimated by a standard one-loop analysis 
Such analysis involves computing the scale q* 
(here q* — 1.446/a), and seeing how the one- loop 
result changes under reasonable variations in q* 
(say 1 < aq* < tt) and in coupling schemes (say 

a Ms(i*) or a v{q*))- 
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Now let us apply what we have learned to the 
calculation of f B in NRQCD. Paralleling Eq. (@), 
the NRQCD expression for fg takes roughly the 
form 

fB<xC 1 (0\O 1 \B) aimQ +—?-(0\O2\B) a , mQ , (3) 

^TOq 

where 

0\ = qjolsQ ; 2 = qjol5p±Q , (4) 

with Q the heavy quark field, q the light quark 
field and ]p ± = 5D i=1 JiDi the spatial Dirac op- 
erator. Note that, unlike in Eq. ([!]), the matrix 
elements in Eq. (^|) depend explicitly on toq . This 
is because toq already appears in the lowest order 
NRQCD Lagrangian. 

As in the HQET case, O2 has a 1/a diver- 
gence proportional to Ox{\ Therefore one expects 
that a renormalon shadow will appear, producing 
large errors in the renormalized 1 /toq matrix ele- 
ment unless the perturbative calculation is taken 
to high order. However, because the matrix ele- 
ment of 0\ also depends on ttiq, the renormalized 
O2 matrix element in NRQCD is not the entire 
0(l/mg) effect. This is a crucial difference from 
the HQET case. The importance of the renor- 
malon shadow here is thus a numerical question. 

The discussion of fg in the Fermilab formalism 
is similar. In Eq. (||), one just replaces 1/(2toq) 
by adi, where d\ is the coefficient of the "rota- 
tion" defined in p6| and is a function of arriQ . As 
arriQ — > 00, d± — > 1/(2citoq), thereby reproducing 
NRQCD. Other the other hand, dx -> as a -> 0. 
This means that the subtraction here does not di- 
verge (and in fact vanishes) as a — ► 0. Note how- 
ever that the relative perturbative uncertainty in 
the renormalized 1/toq operator is unaffected by 
the presence of ad\, which is just an overall fac- 
tor. O2 itself is still power divergent and a renor- 
malon shadow should still appear. Still, the factor 
ad\ may reduce the numerical importance of the 
shadow. 

To discuss the renormalon shadow in fg quan- 
titatively, I start with the Fermilab formalism, 

1 Note, however, that the separation of C\ into c\ and 
ci is now no longer meaningful since the since the loop 
diagrams can generate arbitrary functions of aM . 



using MILC data at = 6.15 with the non- 
perturbative J31J value of the clover coefficient 
Csw- The perturbative corrections have been 
calculated Q at one loop. I define the subtrac- 
tion here as Oi times the difference between its 
complete perturbative coefficient and the coeffi- 
cient when 02 is omitted from the calculation 
p3[] . Since q* has not been computed for this 
quantity, I choose the value that comes from the 
static-light Z A @ (q* = 2.85/a) and use a v (q*) 
fjOfl . To get the uncertainty in the subtraction, 
I then replace q* — > 1/a. (The change is smaller 
when q* is increased, even to 10/a, or when the 
scheme is changed to MS or "boosted perturba- 
tion theory" pq|.) The ratio of the uncertainty 
to the renormalized 1 /toq matrix element is then 
quite large: 60%. I take this as evidence that a 
renormalon shadow is indeed present. 

However, since the renormalized 1/toq opera- 
tor contributes only a small amount to fg, the 
ultimate effect of the the renormalon shadow is 
small. Figure || shows the separate effects of the 
bare 1/rriQ operator and the subtraction, as well 
their combined effects (the renormalized 1/toq 
operator). Clearly, most of fg, as well as its 
rtiQ dependence, comes from the matrix element 
of Ox- The matrix element of the bare 
operator adds 4.6% to fg; while the subtrac- 
tion is —2.6%. The renormalized 1/toq oper- 
ator thus contributes only (2.0 ± 1.2)%, where 
I've included the uncertainty from the renormalon 
shadow. The uncertainty (~ 2 MeV) is smaller 
than many other systematic effects in current 
computations. Note that the error which would 
be made from including the bare 1/toq opera- 
tor without the one-loop subtraction is probably 
larger than leaving out the operator completely. 

In the NRQCD case, most of the numerics (for 
fg B at (3 — 6.0) can be found in Ref. pq |. Here 
the bare 1/toq contributes ~ 11%; the one-loop 
subtraction, ~ —7%. The inclusion of the renor- 
malized 1/toq operator is thus a (4 ± 4)% effect, 
where the error again comes from my varying g* 
in the subtraction. It is no surprise that the ef- 
fect and error here are larger than in the Fermilab 
case: 1/(2q,toq) is considerably larger than d± for 
relevant values of arriQ. However, the 4% un- 
certainty (~ 7 MeV) is still smaller than several 
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Figure 3. 

Effect of bare and renormalized 1 / ttiq operator on 
heavy-light decay constants in the Fermilab for- 
malism. Mhl is the mass of a heavy-light pseu- 
doscalar meson, and Jhl is its decay constant. 
MILC data (/? = 6.15, C S w = 1-644) {u| is used. 



other systematic effects. 

My conclusion is that the calculation of lep- 
tonic decay constants of B mesons is under control 
for both the Fermilab formalism and NRQCD, 
despite the presence of renormalon shadows. The 
issue however needs to be considered on a case- 
by-case basis. There is no guarantee that renor- 
malon shadows are negligible for other physical 
quantities, systems, or orders in perturbation the- 
ory. Indeed, in many computations the perturba- 
tion theory is only put in at (tadpole-improved) 
tree level. The subtractions are therefore left 
out entirely. From the above discussion it seems 
likely that the inclusion of higher dimension op- 
erators in such a situation is worse them omitting 
them entirely, and may lead to large renormalon 
shadow uncertainties (see Sec. ||). 

Two further comments: (1) My numerical es- 
timates of renormalon shadow uncertainties are 
completely standard and do not involve the de- 
tails of the renormalons at all. Therefore it is 
possible to take the term "renormalon shadow" 
as merely a convenient shorthand for the phrase 
"large perturbative uncertainties in power diver- 
gent subtractions." (2) The error I estimate for 
renormalon shadows is of course not the only per- 



turbative uncertainty: many of the terms in C\ 
would be present even if O2 were omitted. Once 
it is ascertained that renormalon shadows are un- 
der control, it is probably preferable to quote 
simply an overall perturbative uncertainty, ob- 
tained by varying q* everywhere. At the current 
state of the art (one loop) the overall perturbative 
uncertainty is generally larger than what I esti- 
mated above for the renormalon shadow. At high 
enough order, my method would probably over- 
estimate the perturbative uncertainty because I 
have not taken into account the expected can- 
cellation of renormalons between c\- and ci-type 
terms. 

2.2. Discretization Errors in MILC Data 

In MILC calculations of quenched fg (|| using 
Wilson and static heavy quarks, the most impor- 
tant source of systematic error was the contin- 
uum extrapolation. This error was manifest in 
the difference between an linear extrapolation of 
f B from all data sets (0.2 GeV _1 <a<0.8 GeV" 1 ) 
and a constant extrapolation from only the finer 
lattices (a<0.5 GeV" 1 ; (3 > 6.0). As of June, 
1999 (with additional running from what ap- 
peared in ||), the former extrapolation gave f B = 
154(11) MeV; while that latter, 180(10) MeV. On 
the basis of the behavior of /bJ/b, we chose the 
linear extrapolation for the central value and took 
the large 26 MeV difference as the discretization 
error. 

Since that time, two developments have signif- 
icantly reduced the discretization error we quote. 
First of all, Tom DeGrand and I have recalcu- 
lated Q the one-loop scale q* appropriate to 
the static-light axial current renormalization con- 
stant Z^ L . Instead of Hernandez and Hill's result 
q* = 2.18/a Hi for tadpole-improved light Wil- 
son fermions, we find a q* that is mildly depen- 
dent on am,Q (as Z^ L is) and is ps 1.4/a for typical 
values of arriQ. We differ from Ref. |37| because: 
(a) we define q* using the complete integrand for 
Z^ L (including the continuum part, which gives 
the toq dependence), and (b) we do not discard 
pieces of the Z^ L integrand which vanish by con- 
tour integration — such pieces do not vanish for 
q* because of the additional \n(q 2 ) factor. Differ- 
ence (b) is responsible for most of the discrepancy. 
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Figure 4. 

Updated MILC data Jlfj] for quenched fs vs. lat- 
tice spacing. Wilson and nonperturbative clover 
( "NP" ) fermions with the Fermilab formalism are 
used. "NP-tad" and "NP-IOY" represent two dif- 
ferent ways of renormalizing the heavy-light axial 
current — see text. 



In the MILC calculation, the static-light q* de- 
scribed above is employed as the central value 
of q* not only in the static-light renormalization 
but also in the (propagating) heavy-light renor- 
malization, for which Za, but not q* , has been 
computed p8[ . The new value q* sa 1.4/a consid- 
erably reduces the a-dependence of Jb with Wil- 
son quarks The difference between the two 
extrapolations of this data is now 15 MeV instead 
of 26 MeV (see Fig. @). 

In addition to reanalyzing the Wilson data, we 
recently completed running at (3 — 6.0 and 6.15 
with nonperturbative |3l|| clover fermions and 
the full Fermilab formalism through order I/ttiq. 
The 1 Iraq operator is always renormalized at one 
loop 



as in Sec. 2.1. For the rest of the 



renormalization of the heavy- light axial current, 
we use either the one- loop calculation |3^j , or a 
"tadpole renormalization" Ansatz designed to re- 
produce the nonperturbative renormalization of 
Rcf. Bit] at small mass and have a sensible limit 



as itiq — > oo. These two approaches are shown as 
"NP-IOY" and "NP-tad," respectively, in Figji. 



NP-tad is described in more detail in Sec. 2.3. 
Both NP-IOY and NP-tad are correct to 0(a). 
Higher order effects (e.g., 0(a 2 ), 0(g 4 )) are com- 
plicated but, hopefully, rather small. Each ap- 
proach is fitted to a constant; the difference in 
extrapolated value from the two NP approaches 
gives an estimate of these higher order effects. 

From Fig. [| it now seems clear that it was a 
mistake to choose only the linear extrapolation to 
find the central value of the Wilson data in the 
continuum: the linear fit gives the lowest extrap- 
olated value of all four approaches. (In the case 
of /b s //b, the linear extrapolation differs even 
more from the other three fits.) MILC currently 
averages all four approaches for the central value 
and defines the discretization error as the stan- 
dard deviation of the four. The result appears 
as the quenched "MILCOO" point in Fig. [I] and 
Table 0. 

2.3. "Nonperturbative" Heavy-Lights 

APE Ull' UK Q CD ©@: and MILC © 

have simulated heavy-light physics using nonper- 
turbative clover fermions. The starting point is 
the expression for the improved, renormalized ax- 
ial current A$: 



Aq =Z A ^An q K Q (l 
Aq = 97o75 Q 5 P5 



bAamQ$ 



2 

975 Q 



)[A +c A adoP 5 



(5) 



where K q , kq are the light and heavy quark 
hopping parameters, respectively, and uiq^q is 
the bare heavy quark mass (arnQ^ = 1/(2k) — 
1/(2k c )). For simplicity, the light quark mass has 
been set to zero. 

All three groups use the nonperturbative val- 
ues of Z A and c A in Eq. (||) (as well as the clover 
coefficient Csw) computed by the Alpha collab- 
oration The groups differ, however, on the 
choice of the coefficient 6^, which was not com- 
puted nonperturbatively in Ref. pl| . At (3 = 6.2, 
for example, APE uses bA = 1-24 from the one- 
loop calculation and boosted perturbation 
theory @; while UKQCD uses b A = 1.47 from a 
preliminary nonperturbative calculation of Bhat- 
tacharya et al. |39|. This difference accounts for 
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any disagreement between APE and UKQCD re- 
sults that remains once the scales are set in the 
same way. MILC's bA is taken from perturbation 
theory |40|j , but with coupling ay[q*), with q* 
chosen as the value (ss 1/q) which produces the 
nonperturbative result |31| for the similar quan- 
tity by- This gives values of 6a (1.47 at [3 = 6.15 
and 1.42 at /3 = 6.0) that are quite close to those 
of UKQCD .0 

APE and UKQCD apply Eq. © directly for 
moderate values of am,Q$ (up to ~0.5 at = 6.2 
and ~ 0.75 at j3 = 6.0). This allows them to 
reach a maximum meson (pole) mass of ~ 2 GeV 
at (3 = 6.2, where I've set the scale by m p . Using 
HQET, which implies that /hlVMhl should be 
a polynomial in 1/Mhl (up to logs), they then at- 
tempt to extrapolate the results up to the B mass. 
(Mhl and Jhl are the mass and decay constant 
of a generic heavy-light pseudoscalar meson.) 

There are of course systematic errors associ- 
ated with this approach. First of all, one does 
not know a priori what order polynomial to use; 
there is a large difference (~ 10%) between the 
results of linear and quadratic extrapolations. A 
second systematic effect is more subtle. Although 
discretization errors in Eq. (|^) are very small for 
o-rriQfi <C 1 and appear to remain quite small 
even for the maximum arriQfl used by ULQCD 
and APE, those errors grow rapidly with arriQfi. 
Indeed, for arriQ o — + oo, Ar in Eq. (||) does not 
approach a static limit as it should, but goes to 
— oo because ca < and do-Ps ~ sinh(M.Hx,i) ~ 
o- m Q,o- (Here Mhl,i is the meson pole mass.) 
Even if ca were zero, Ar would still blow up 
because of the term bACimQfi. Thus, small dis- 
cretization errors may be magnified significantly 
by the extrapolation to the B. 

To estimate the latter error, one may compare 
to the MILC "NP-tad" approach g§. The MILC 
goal is to replace Eq. (||) by an expression that 
is equivalent through O(a) but which gives a sen- 
sible limit for fn l V Mhl as arriQo — > oo. The 
result is then used at arbitrary amQ^, a la Fer- 



milab. To do so, we first define 
(0\d P 5 \HL) 
mQ,o{0\A o \HL) 

where HL is a generic heavy-light pseudoscalar 
meson. Due to a cancellation of sinh(MHx,i) 
(from do) and the explicit toq,o in the denomina- 
tor, one expects R has a finite limit as amg^o — ► 
oo. This is confirmed by simulations. Then 



2 At j3 = 6.0 and 6.2, a new nonperturbative calculation of 
all the coefficients, including bA, is now available ]4l| . It 
will be interesting to see how these values affect the results 
for f B . 



A£ = Z A y/iK q K Q ^l + (b A + 2c A R)am Q:0 A (7) 

gives results for (0\Aq\HL) that are identical 
to Eq. (||) through 0(a). However, because 
KQaniQfl — > 1/2 as kq — > 0, Eq. (Q) has a static 
limit, unlike (||). Indeed, (^) is just a version of 
the Fermilab formalism at tadpole-improved tree 
level (the 1/toq operator is omitted for simplic- 
ity), but with a special (mass-dependent) value 
for the tadpole factor: uq = (bA + 2caR)~ 1 - 
The similarity to tadpole improvement (within 
the context of nonperturbative renormalization) 
is the reason for the name "NP-tad." 

Figure || shows the effect of reanalyzing the 
UKQCD results with the NP-tad normalization. 
This O(o?) effect is 2% to 5.5% on their data 
points, but 8.5% after extrapolation to the B. 

An additional difference between MILC and 
UKQCD or APE is that MILC, following the 
Fermilab formalism, defines the meson mass as 
the kinetic mass M2, rather than the pole mass 
Mi. Since the difference in masses is 0(a 2 ) this 
does not affect the O(a) improvement. M2 is de- 
termined by adjusting the measured meson pole 
mass upward by AM = mi — mi , where m2 and 
mi are the heavy quark kinetic and pole masses, 
respectively, as calculated in tadpole-improved 
tree approximation pq , with the mean link fixed 
by k c . Figure [5] shows the effect of this shift on 
the UKQCD data (after first changing to NP-tad 
normalization). One can check the MILC deter- 
mination of M 2 by using instead the M 2 values di- 
rectly computed by UKQCD from the meson dis- 
persion relation. The resulting curve is not shown 
in Fig. ^ because it is indistinguishable from the 
fit to the diamonds. 

Note that the effects of the NP-tad norm and 
the shift by AM almost cancel, so that the fi- 
nal results of UKQCD and MILC are actually 
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Figure 5. Effects, at f3 = 6.0, of various 
treatments of heavy-lights using nonperturbative 
clover fermions. UKQCD results (crosses) have 
been adjusted to same scale (from /„■) as MILC 
results (diamonds). Squares show the effect of 
adjusting the crosses using the NP-tad normal- 
ization. Octagons include the tadpole-improved 
mass shift (described in text) on the squares. 



quite consistent as long as the scale is set in the 
same way. Since the two effects are logically inde- 
pendent, however, the 17 MeV difference for Jb 
among the curves in Fig. |5| is a measure of the 
discretization error. This error is still 15 MeV at 
(3 = 6.2. I would not, however, argue that this er- 
ror should be added to those quoted by UKQCD. 
They already include a discretization error error 
obtained by comparing the results at (3 — 6.2 and 
6.0. Indeed, the discretization error computed 
that way is 16 MeV (using an m p scale), identical 
to the estimate above. 

The goal of the approaches discussed in this 
subsection is to treat the heavy-lights in a way 
that takes advantage of the known nonperturba- 
tive renormalizations for light-lights. But since 
one cither works in or extrapolates to a region 
where arriQ is not small, the systematic errors are 
not necessarily small. Furthermore, because the 
methods agree at 0(a), there is no way a priori to 
distiguish among them at the level of Eq. (|J) . If 
one insists on using the nonperturbative informa- 



tion for the heavy-lights, then I prefer the NP-tad 
approach (and equating the physical mass with 
M 2 ) because (a) having a static limit seems de- 
sirable if we want to use (or extrapolate to) large 
arriQ, and (b) it appears to scale somewhat bet- 
ter: the difference between Jb at (3 = 6.0 and 6.15 
is only 2 MeV (see Fig. ||). However, this is far 
from definitive: point (a) is subjective, and point 
(b) is not very convincing with only two lattice 
spacings. 

The situation will improve soon. Although a 
true nonperturbative treatment of heavy-lights 
seems difficult, a nonperturbative computation 
of the 0(a) renormalization/improvement coef- 
ficients for the static-light axial current is close 
to completion p3[ |. This will be an important 
advance, since it will allow one to compute fs 
with an interpolation between two nonperturba- 
tive calculations. The extrapolation error of the 
UKQCD or APE approach will therefore be much 
reduced. One will also be able to see how well 
NP-tad is really doing in the intermediate region. 

A alternative computation with nonperturba- 
tive clover fermions that can be done today is 
the "NP-IOY" approach mentioned in the previ- 
ous section. However, since it is just a marriage 
of standard techniques, it is perhaps misleading 
to give it a special name: One simply confines 
the use of the nonperturbative information to the 
light-light sector, where it is completely justified. 
In this case that means setting the scale with, for 
example, or m p . The straightforward Fermi- 
lab approach, with one loop perturbation theory 
p2|] , is then used for the heavy-lights. 

It is also worth mentioning here the calcula- 
tion of Ref. @. This is similar to Refs. (l|Jl|] in 
that the B meson is reached by extrapolation from 
relatively small values of arriQ. However the im- 
provement is done perturbatively, and the EKM 
normalization |2(| (but not the shift Mi — > M-i) is 
used for the central values. I believe that the er- 
rors quoted in Ref. adequately account for all 
systematic effects (within the quenched approxi- 
mation) of that calculation. 

2.4. Unquenching fs 

Figure ^ shows the current status of fs cal- 
culations with two flavors (Nf = 2) of dynam- 
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Figure 6. 

World data for fs with Nf = 2 vs. lattice spac- 
ing. Only statistical errors are shown. "MILC99" 
is Ref. see Table [j]for other references. 



ical (sea) fermions. One improvement over last 
year is that the very low value from the MILC 
"fat-link" clover computation has now been 
understood. From measurements of the quenched 
static quark potential, we have shown that fatten- 
ing smooths out the potential well at short dis- 
tances relative to the corresponding "thin-link" 
(standard) case. Indeed, for the fattening MILC 
used (10 iterations of APE smearing with sta- 
ple coefficient a = 0.45), the potential at distance 
1 (for f3 = 5.85) is increased by 45%. At distance 
\/3, the increase is only to 8%. And by distance 
3, it is under 1%. Such an effect is not surprising; 
the expected range for this amount of fattening 
is ~ y/10 x 4.5 pa 2.1 lattice spacings (j46j. One 
expects that when the potential well at short dis- 
tance becomes shallower, decreases, since it is 
basically the wave function at the origin. 

On quenched configurations at (3 = 6.0, 
/bV/b' 11 = 0.76 when "thin" is the NP-tad ap- 
proach, and 0.73 when "thin" is NP-IOY. The 
extent of the effect was somewhat surprising to 
us, since the light hadron spectrum behaves quite 
well with fat-link clover fermions jlTj and since 
the potential is only changed by an amount > 10% 
for distances < 0.16 fm. The difference between 



fat and thin may be exaggerated by the fact that 
in the fat case the one-loop renormalization of 
heavy-light Za has not been computed and the 
static-light Za [|6|3j] was used instead. Despite 
this, the ratio /Jj at // B hin should be roughly the 
same at Nt = 2, am = 0.01, (3 = 5.6 (the pa- 
rameters MILC used in the unquenched case) as 
at quenched = 6.0 because the lattice spac- 
ings and couplings are quite close. Under this 



/•fat, corrected 
lB,N f =2 



= f. 



fat 
B.N, 



f Blench/ f Blench * Fig. | with thin=NP-tad. 



assumption, I plot f E 

p thin / J'fat 

B, quench/ J B,q\ 

The MILC result for Jb,n 1= 2 in the continuum 
is obtained by extrapolating the Wilson data un- 
der two different assumptions that seem likely to 
bracket the real behavior: (1) constant behavior 
as a — > 0, and (2) a linear decrease from small- 
est current a values («0.45 GeV -1 ) with a slope 
equal to that of the linear fit to the quenched 
data (see Fig. [|). The average of the two proce- 
dures is taken as the central value; the standard 
deviation, as the extrapolation error. The (pre- 
liminary) result shown in Fig. |l| and Table |l| then 
is: f B = 191(6)(lig)(+ ( l ) 1 ) MeV, where the errors 
are respectively statistical, systematic within the 
Nf — 2 approximation, and an estimate of the er- 
ror of neglecting the strange sea quark. Although 
the consistency of the corrected fat-link results 
with the Wilson data is comforting, the former is 
not included in MILC final results because there 
is no guarantee that the correction factor is the 
same in the quenched and full cases. 

As in the quenched case, the preliminary CP- 
PACS results for fg with Nf = 2 use the Iwasaki 
gauge action ||. They find f B ,N } =2 = 

208(10)(11) MeV with the Fermilab formalism 
(heavy clover) and 204(8) (15) (+44) MeV with 
NRQCD; these results are shown in Fig. The 
smallest-a points in Fig. |] are taken as the 
"continuum" values, and the discretization er- 
ror comes theoretical estimates: 0(a 2 A 2 ) and 
O(aaA) in the heavy-clover case. However, 
my comment about discretization errors in the 
quenched CP-PACS data applies here too. In- 
deed, comparing Fig. || with Fig. ||, one sees that 
the lattice spacing dependence of the Nf = 2 data 
is if anything greater than in the quenched case. 
The theoretical error estimates do not account for 
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the full amount of the variation with a. I there- 
fore think the procedure for finding central values 
and errors is overly optimistic. 

My own rough estimate from a const + a 2 fit 
to the CP-PACS heavy-clover data would be 
fB,N f =2 = 190(12) (26) MeV. The systematic er- 
ror in my estimate is dominated by the discretiza- 
tion error, taken from the difference of the extrap- 
olation and the smallest-a point. Although one 
can certainly argue with my extrapolation (there 
are for example aaA errors in addition to a 2 er- 
rors) I believe the systematic error of ~ 26 MeV 
is reasonable. 

One difference between CP-PACS and MILC 
data shown in Fig. ^| is that the MILC data is 
"partially quenched" — the sea quark mass is 
held fixed while the valence mass is extrapolated 
to the physical u,d mass. The CP-PACS points 
are fully unquenched, with the valence and sea 
masses extrapolated together. However, MILC 
has repeated the analysis in a fully unquenched 
manner; the points are then still roughly constant 
in a, and the final result is changed only slightly. 
The difference is included in the MILC systematic 
error. 

In the NRQCD N f = 2 points in Fig. | 
the systematic errors are quite large. The dom- 
inant source of these errors is the ~ 25% differ- 
ence in the lattice scale from that set by light 
quark physics (m p , say) and the IP- IS T split- 
ting. In previous, quenched simulations, it was 
generally assumed that the scale difference was 
due to quenching and would go away with dy- 
namical quarks. Although the difference is in 
fact somewhat smaller for N 



f 



2 than in the 



quenched approximation, it is still large; there is 
no indication that it would vanish in the physical 
case, i.e., 3 light dynamical quarks. Perhaps the 
errors for heavy-heavy spectra in NRQCD are un- 
derestimated (see Sec. ||). I do not, however, have 
the data and calculations to subject this suspicion 
to a quantitative test. 

Looking at Fig. [l] and recalling my remarks on 
discretization effects, one sees that the system- 
atic errors on in both the quenched and the 2- 
dynamical-flavor cases are still rather large over- 
all. However, comparing quenched and Nf = 2 
results by the same groups, for which many of 



the systematic effects cancel, one can say with 
some confidence that fB,N f =2 is about 10-15% 
larger than fs, quench- The effect of unquenching 
is roughly the same for fg s ; while the increase 
for fr> and /n appears somewhat smaller, about 
3-8% Qll-pl. 



Taking into account all above remarks, I arrive 
at the following world averages: 



,f B , quench 

= 175(20) MeV, 
^) = 1-15(4), 

JB J quench 

f D , = 255(30) MeV, 



f B = 200(30) MeV, 
^ = 1-16(4), 

IB 

(8) 



where quantities without the "quench" subscript 
are supposed to be full QCD quantities, including 
an estimate of the effect (and error) that would re- 
sult from including the strange sea quark. Start- 
ing with the quenched results and adjusting them 
by the expected increase in the full theory would 
result in similar unquenched values and errors. 

3. B PARAMETERS 

Figure fj] presents most of the recent world data 
on the B-B mixing parameter Bb a , where the 
"hat" indicates the renormalization group invari- 
ant quantity at NLO. Three types of calculations 
are shown: heavy-light with relativistic heavy 
quarks (here, clover quarks) extrapolated to the 
B meson mass [ppTI PI], static-light (4§|4|], and 
NRQCD-light gff. 

In the relativistic heavy-light case, the extrap- 
olation to the B is motivated by the heavy quark 
effective theory. HQET implies that Bp d is a 
power series (up to logs) in 1/Mhl, where Bp d 
is the (left-left) B parameter of a heavy-light me- 
son P (mass Mhl) with a physical d quark and 
an arbitrary heavy quark. Both APE [^2j and 
Lellouch and Lin |L9|] make linear fits in 1/Mhl, 
resulting in the extrapolated points shown in 
Fig. 0. Although it must be admitted that the 
data of Ref. jl9| fit the linear form quite well, I 
have concerns about trusting the HQET down to 
M HL ~ l-2GeV (1/M HL w 0.85 GeV^ 1 ), which 
is the lightest mass shown in Fi g. M. I ndeed, since 
the relativistic points in Refs. [ |12||1S| ] suffer from 
some of the same systematic effects discussed in 
Sec. 2.3 {e.g., the choice of Mi or M2), the ex- 
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Figure 7. 

Selection of recent data for quenched Bp d vs 
1/Mhl, where P is a general heavy- light pseu- 
doscalar and Mhl is its mass. Only statistical er- 
rors are shown for heavy- light data and extrapola- 
tions ||@|lj. Static-light points g 
perturbative systematic errors 
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uses NRQCD-light; the 4 points at each mass rep- 
resent different ways of treating the 0(a 2 s ) and 
O(a s p/mo) corrections. 



tremely linear behavior may well be an accident. 
I therefore will omit the 0.85 GeV -1 point from 
my analysis. 

References |l^|li| do not take into account the 
static-light results {e.g., Refs. [ p8p9| ]) in their ex- 
trapolations. APE justifies this by noting that 
the perturbative error on the static-light B pa- 
rameter is large, while the heavy-light points in 
12 1 use nonperturbative normalization. In my 
view, however, the arriQ ~ 1 systematic errors 
in the heavy-light case (sec. 2J5) are compara- 
ble to the perturbative errors of the static-light 
method. Therefore, I choose to make a quadratic 
fit (in X/Mhl) to the APE points as well as the 
Gimenez and Reyes/UKQCD static-light point. 
(The latter probably has the smallest systematic 
errors of the available static-light data, since it 
uses improved light quarks at the weakest cou- 
pling, (5 = 6.2.) My fit is shown in Fig. At 



the B mass, it gives B Bd — 1.30(6) (statistical 
error only). Note that the extension of the fit 
comes reasonably close to the 0.85 GeV -1 point; 
alternatively, including this point in the fit would 
produce only a 1% change in B Bd . 

An important advance in the past year is the 
calculation |5l[| of the one-loop renormalization 
of the AB = 2 4-quark operators in the NRQCD- 
light case. Previously p2[ , the static-light renor- 
malization constants had to be used instead, and 
the result for B Bd was ~20% low compared to 
the extrapolated values of Refs. jl2],[l9) and still 
~ 15% below my interpolated value. (Note, how- 
ever, that the perturbative error in Ref. |5^] was 
estimated rather realistically as ~ 10%.) The 
current perturbative calculation includes terms 
of 0(a s /(amg)), but not those of 0(a s p/mQ), 
where p is the typical quark 3-momentum in the 
B meson. The various squares at the same values 
of 1 /Mhl in Fig. [?| represent an attempt to esti- 
mate the systematic errors of 0[a s pjmQ) as well 
as 0{o? s ). The results are consistent with those 
of my interpolation. Unfortunately, I do not have 
the data to perform a renormalon shadow analy- 
sis on the NRQCD-light results, but my guess is 
that the effect is no larger than in the f B case. 

Taking into account my interpolation and the 
results of Refs. J5^j5l[], as well as my esti- 
mate of the systematic errors in both, I quote 
BB d ,quench = 1.30(12). It is not clear what 
quenching error should be assigned here. From 
quenched chiral perturbation theory (QChPT), 
Sharpe |53| estimates a 10% error. There are 
not many simulations that address this question. 
Gimenez and Reyes [Q find a ~ 11% difference 
between quenched and unquenched results in the 
static-light case. However, the systematic effects 
in the two results are rather different, so it is dif- 
ficult to draw a strong conclusion. MILC also has 
some quite preliminary results p5[ ] in the static- 
light case, suggesting a smaller quenching error, 
but again the systematics in the comparison are 
not under good control. At this point therefore I 
quote the QChPT 10% as the quenching error on 
the "world average" B Bd : 



B Bd = 1.30(12) (13); f B J B Bd =230(40) MeV, (9) 
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where the uncertainty on fg y Bs d is supposed to 
include all sources of error. Since QChPT for fs 
]5C| seems to overestimate somewhat the quench- 
ing errors, I regard the assumed 10% quenching 
effect on Bb a as quite conservative. 

For the ratio Bb 3 / Bs d , all groups get a number 
very close to 1. Furthermore, the experience with 
Jb 3 //bj as well as a preliminary simulation (55) , 
suggest that the quenching errors on the ratio are 
considerably smaller than the ~ 5% inferred from 
QChPT ||§. My world averages are then: 



|^ = 1.00(4); 

Bn, 



1.16(5), (10) 



f Bd \B B 



where I have used Eq. (g). All sources of error, 
including quenching, are supposed to be repre- 
sented in Eq. @. 

Before concluding this section, I wish to men- 
tion results for B$, the B parameter for the 
scalar-pseudoscalar AB — 2 operator, which is 
relevant for the Standard Model prediction of the 
width difference AT/T for B s mesons. B$ has 
been computed with heavy quarks that are rel- 
ativisticjl^HIJS nonrelativistic |^j5|] and 
static ]54 |. The last reference includes a first 
look at unquenching effects. The results from all 
the groups appear consistent with Bg IS (m.i,) = 
0.82(8), where I have defined B~fs as in (5§, and 
have tried to include all sources of error. There 
are however large 1 /rriQ corrections as well as dis- 
agreements in how to go from Bs to the width dif- 
ference, so more work is required before we have 
a Standard Model prediction for AT bJT b,- 



4. SEMILEPTONIC FORM FACTORS 

New work this year on semileptonic decays has 
dealt almost exclusively with B — > irii/, so I will 
focus on that process only. The hadronic matrix 
element can be parameterized with form factors 
/+ and /o defined by 
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Figure 8. 

Recent data for B — > -nlv formfactors /+ and fo- 
Only statistical errors are shown. 



where Vn is the vector current and q = p — k is 
the 4-momentum transfer to the leptons. 

Figure || shows recent results for the form 
factors. Various approaches have been triedj^] 
UKQCD ||,|l| and APE @ use relativistic 
quarks and extrapolate to the the B; JLQCD |(i(J 
treats the heavy quark with NRQCD; and the 
Fermilab group J6l| uses their formalism to sim- 
ulate with the heavy quark at the b mass. 

In all the approaches, the first step in analysis 
of the raw lattice data for the form factors (or 
for (ir(k)\V fl \B(p)) in the Fermilab case) is an in- 
terpolation or fit. For APE or UKQCD, this is 
a fit of /_)_ and /o to expected phcnomcnological 
forms. APE uses the nice "B-K" form j6j|, which 
automatically incorporates the known scaling be- 
havior in I/Mjil at q 2 = q^ lax (the end point) 
and at q 2 = 0, and enforces the kinematic con- 
straint /+(0) = /o(0). UKQCD employs more 
conventional forms such as pole/dipole, although 
fits to the B-K expressions do appear later in their 
analysis. In practice, what is important at this 
stage is that the fit smooths out data which is 
inherently rather noisy. For example, I refer the 
reader to Figure 2 in Ref. [0, which I believe is 
typical of the kind of data seen by all the groups. 



(11) 



3 A form factor study using coarse anisotropic lattices |62J] 
is not far enough along yet to be included in this compar- 
ison. 
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One must thus keep in mind, when looking at 
plots like Fig. ||, that the apparent smoothness 
of the data is a result of a procedure that starts 
with an "averaging" process over a range in q 2 
and therefore produces highly correlated points. 

The next step is a chiral extrapolation. In 
the heavy meson rest frame, q 2 — Mjj L + M 2 — 

2Mjjl \J M 2 + k 2 . This means that an extrap- 
olation of form factors at fixed, small k should 
contain a term linear in M Vl through the implicit 
q 2 dependence, in addition to the usual powers 
of M 2 Q. Alternatively, one can perform the 
extrapolation at fixed q 2 |5£j or fixed v ■ k pl| , |6C| 
(v is the heavy meson 4- velocity), or simply ex- 
trapolate the B-K parameters (or similar ones) 
rather than the form factors themselves (2lJ . Yet 
another choice is to avoid the small k region al- 
together (the decay rate is after all small there) 
and make a standard chiral extrapolation at fixed 

UKQCD and APE then need to extrapolate 
in 1/Mhl to the B mesons. UKQCD extrapo- 
lates /+ and /o at fixed v ■ k, assuming that q 2 is 
close enough to q 2 lax that the HQET forms valid 
near the end point may be used. APE extrap- 
olates the B-K parameters themselves (and also 
repeats the UKQCD-type analysis as a check). 
The two extrapolations are consistent within er- 
rors (see Fig. ||) , suggesting that the HQET end- 
point forms are indeed appropriate. Of course, 
just as for fs and Bb, either extrapolation in- 
troduces two sources of systematic error: (1) un- 
certainty in what order polynomial (in 1/Mbl) 
one should use, and (2) magnification of small 
0(a 2 Mft L ) errors. Both groups do an adequate 
job of estimating (1). However, I think that error 
(2) could use more study along the lines of Sec. 


JLQCD and the Fermilab group work in the 
region of the B mesons, so extrapolation in the 
heavy mass is not required. Although their re- 
sults in Fig. U are similar, there may be some 
quantitative disagreement, especially in the nor- 
malization of /o and in the q 2 dependences of both 
form factors. I have some concern that the Fermi- 
lab data for / + does not appear to show the ex- 
pected pole-like rise toward the end point seen by 



the other groups. (Remember that each group's 
points are highly correlated, so that the shape of 
the q 2 dependence is probably highly significant.) 
Of course, for APE and UKQCD, this behavior is 
put in ab initio by the fit to the phenomenolog- 
ical form. JLQCD, on the other hand, does not 
appear to be forcing this behavior by their anal- 
ysis procedure, so the fact that it does come out 
seems encouraging. 

Since the Fermilab and JLQCD (NRQCD) ap- 
proaches are fundamentally quite similar in the B 
meson region, it is not easy to guess the source 
of the differences. One possibility is discretiza- 
tion error. The Fermilab group (alone among the 
four) extrapolates to the continuum from a range 
of a (P = 5.7, 5.9, 6.1), while JLQCD works only 
at (3 = 5.9. However, the a-dependence seen by 
the Fermilab group is mild, and it seems unlikely 
to me that this is the explanation. 

Perhaps a more likely culprit is the current 
renormalization. JLQCD is working to one-loop 
order in perturbation theory. The Fermilab group 
does much of the renormalization nonperturba- 
tively (by comparing to the known normalization 
of the diagonal vector current), but the remain- 
der is still put in at tadpole-improved tree level 
for the data in Fig. |[ They will be including the 
complete one-loop calculation shortly, and they 
assure me that it makes little numerical differ- 
ence. It will, however, be interesting to investi- 
gate the size of the renormalon shadows for both 
the JLQCD and the Fermilab computations. 

Another issue that has bedeviled lattice B — > 
■niv calculations for several years is the soft-pion 
relation (SPR) at the end point, foiVmax) = 
/b//tt p5| , |66[ . Some early calculations found vi- 
olations by as much as a factor of 2 |67]], and 
although the situation has improved somewhat, 
the problem has not gone away. 

At first sight, it seems surprising that the end- 
point region should give trouble. After all, lattice 
calculations are easiest at the end point where no 
3-momentum insertion is required (k — in the 
B rest frame). However, the SPR is only guar- 
anteed to be valid in the limit M„ — > 0. As dis- 
cussed above, the extrapolation in M v when k is 
fixed and small is tricky because q 2 changes and 
has a term linear in M r oc y/m quar k- APE and 
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UKQCD get consistency with the SPR using chi- 
ral extrapolations that avoid sitting at the end 
point (instead fixing q 2 or v ■ k). The end point 
is then reapproached after extrapolation to the 
physical q^ ax and the B meson mass. But, after 
all the extrapolations, the systematic errors are 
rather large, so the agreement with the SPR is 
not very convincing. Indeed, in the case of APE, 
the SPR is used as a criterion to choose whether 
one should do linear or quadratic fits in 1/Mhl, 
so the agreement cannot then be seen as a success 
of the computation. 

JLQCD, on the other hand, studies the SPR 
in detail by sitting at the end point with high 
statistics (2150 configurations). Figure || shows 
the chiral extrapolation of fo{q 2 lax ). Note that 
the fit (solid line) which takes into account the 
implicit linear dependence on M n comes closer to 
/b//tt than the purely quadratic fit in M n but 
still misses by a wide margin. This illustrates 
how hard it will be to verify the soft-pion rela- 
tion directly: the linear dependence on M T causes 
a sharp rise, but only at very small mass where 
there is no data. This behavior is expected from 
the effects of q 2 poles in the form factor in the 
unphysical region just beyond the end point. A 
global fit to the functional dependence on v ■ k and 
Mn expected from the HQET-suggested form fac- 
tors of Burdman et al. |(3(| , does somewhat bet- 
ter, but still does not give convincing agreement 
with Jul f-n- The Fermilab group has indepen- 
dently decided to use the form factors of Ref. |6(| 
to study the end-point behavior, and it will be 
very interesting to see how their results compare 
to those of JLQCD. 

Of course, one can always take the attitude that 
the SPR is irrelevant to phenomenology because 
there is no rate near the end point and, besides, 
it is /+, not /o, that determines the rate into ir/iv 
or Trev. However, a credible demonstration of the 
SPR would be very helpful in convincing the non- 
lattice community that we have the form factor 
computation under control. 

5. SPECTRA 

Although I do not have space for a detailed 
review of recent spectral calculations for hadrons 



V-4 



A global fit 



O (am ) + (am ) 
□ (am f 



0.02 0.04 

a m 



Figure 9. 

Extrapolation of /o at the end point to the chiral 
limit by JLQCD [ p0| . The open triangle comes 
from a global fit to the form in Ref. fifl . 



with one or more heavy quarks (Refs. |68|-p0| ) , I 
will make some remarks about a few features I 
find interesting and/or confusing: 

(1) It seems clear that the computed quenched 
hyperfine splittings for heavy-light mesons are 
considerably smaller than the experimental val- 
ues. For example, Mb* — Mb 3 in the quenched 
approximation appears to be about half its phys- 
ical value |38 71 7%, although some older calcu- 
lations e.g., Ref |B1|, found smaller (~ 20%) dis- 
crepancies. For the D s system, the effect is also 
present at about the 20% level in most Jn] , |32"|]8l| ] 
calculations (with one recent exception]70(|). The 
sign and, very roughly, the magnitude of the ob- 
served discrepancy is expected from arguments 
about the effects of quark loops on the short dis- 
tance potential. However, I know of no direct 
evidence from simulations that quenching is the 
culprit. 

(2) There is a disagreement between Refs. J7^] 
and |6Sl l in the computed P-wave splittings for B 
mesons (Mb* — Mb* or Mb* — Mb, or the cor- 
responding quantities for B s ). The disagreement 
is considerably larger than the systematic errors 
quoted by both groups and deserves further study. 
My guess is that the problem lies in excited state 
contamination: These mass differences are quite 
noisy, and it is seems very difficult to identify a 
"true plateau" in the effective mass plots. 
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(3) For NRQCD computations of splittings in 
charmoniurn, Stewart and Koniuk [fTSt confirm 
previous evidence Q that the results are very 
sensitive to the order (in the velocity v) to which 
the Hamiltonian is corrected and how the tad- 
pole improvement is done (plaquette or Landau 
link). Indeed, including the 0(v 6 ) terms, with the 
renormalization done at tadpole-improved tree 
level, moves the computed hyperfinc splitting sig- 
nificantly further from experiment. I believe this 
is a "smoking gun" for renormalon shadow ef- 
fects: At tadpole-improved tree level the presence 
of higher dimensional operators has no effect on 
the coefficients of lower dimensional operators — 
there is no subtraction of power law divergences. 
Just as for /b, putting in the higher dimensional 
operators without the subtraction is then worse 
than leaving out those operators entirely. Unlike 
the Jb case, however, the effect here appears to 
be numerically dramatic. 

(4) In the bb system, NRQCD is clearly better 
behaved than in cc, since the relativistic correc- 
tions are smaller. It is therefore possible to find 
convincing evidence for quenching effects ||77| , [73| 
if one is careful to keep other systematic effects 
(lattice spacing, order in v, type of tadpole im- 
provement) fixed when comparing quenched and 
unquenched simulations. However, even for bb, I 
believe it is important to reduce renormalon shad- 
ows by doing at least one-loop renormalization 
of the Hamiltonian before one can make reliable 
comparisons with experiment for most quantities. 

(5) Following pioneering work by Klassen [ pi[ , 
there have been several recent simulations of char- 
moniurn on quenched "anisotropic relativistic" 
lattices j76|-[78[. The anisotropy £ is defined by 
£ = a s /a t , where a s and at are the spatial and 
temporal lattice spacings, respectively. Choos- 
ing £ > 1 (typically 2 < £ < 5) allows one to 
keep atvrtQ <C 1 and avoid large mass dependence 
in the improvement coefficients. Effectively, the 
heavy quark is relativistic.^ Although the im- 
provement coefficients through 0(a) can in prin- 
ciple be evaluated nonperturbatively, in practice 
the spatial and temporal clover coefficients are 




^Experiment 

• CP-PACS, ^-3, plaquette, r 
■ CP-PACS, i=3, plaquette, 1P-1S 
O Klassen, ^=3, Landau, r D 
AChen, ^=2, Landau, r 

♦ CP-PACS, ^-3, Landau, r 



0.020 0.040 0.060 

a, 2 (fm 2 ) 



Figure 10. 

Hyperfine splitting in charmoniurn vs. a 2 s for 
relativistic anisotropic lattices. Data is from 
Refs. figure is taken from @. 

determined at tadpole-improved tree level, either 
by the Landau link or the plaquette. 

Figure [l^ shows the continuum extrapolation 
for the charmoniurn hyperfine splitting for various 
choices of anisotropy, tadpole factor, and scale 
determination. The good news is that the results 
from this method seem quite precise compared 
to previous relativistic calculations on isotropic 
lattices and do not have the sensitive de- 
pendence on the tadpole factor seen (7|j83|] in 
NRQCD. The bad news is that there appears to 
be some remaining dependence on the tadpole 
factor in the continuum limit (compare the open 
triangles and filled circles).^ The reason for the 
discrepancy is not yet clear. Additional data at 
larger a s with the plaquette tadpole may be use- 
ful in sorting out what is happening. 

(6) The anisotropic relativistic approach has 
also been quite helpful in getting good signals for 
quarkonium hybrids [f78| . As in glueball calcula- 
tions , anisotropy gives more usable time slices 
before the signal is lost in noise. Another recent 
hybrid calculation |]9| also employs anisotropic 
lattices, but in the NRQCD context. 

(7) Kronfeld |7£| has recently argued that a 



4 To see that a a M ~ 1 does not introduce uncontrolled 
maps dependence requires an analysis along the lines 



5 Differences in the continuum value due to different scale 
choices (ro or the IP-IS splitting) are presumably ex- 
plained by quenching and/or phenomenological uncertain- 
ties in ro (see Sec. CJ). 
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heavy quark on the lattice can be represented by a 
continuum HQET where the dependence on arriQ 
is completely absorbed into the short-distance co- 
efficients of the heavy quark operators. The oper- 
ator matrix elements then only have mild lattice 
spacing dependence though the light degrees of 
freedom. This nice insight is then applied Q to 
the heavy-light meson spectrum on the lattice in 
an attempt to compute the HQET parameters A 
and Ai. (A is the meson binding energy in the 
static limit; — Ai/2mQ, the heavy quark kinetic 
energy.) 

I have some concerns about Ref. |80 . One 
could find the continuum HQET by a two-step 
process: (a) for fixed a in the lattice theory, take 
rriQ large and arrive at a lattice HQET, and (b) 
order by order in perturbation theory, replace the 
lattice regularized HQET operators by their con- 
tinuum counterparts. Although this not the pro- 
cedure employed in [Q, I believe the approaches 
should be equivalent. If I am right, then there 
are effective power law divergences and renor- 
malon shadow effects^ in the coefficients of the 
continuum operators (resulting from step (b)). 
Therefore, a one-loop calculation may not be suf- 
ficiently accurate to extract A and Ai. Indeed, 
parameterizing the discretization errors on A as 
A(a) = A(0)(1 — aM), the slope M. has a value 
of sa 0.65 GeV. Yet, if the lattice spacing de- 
pendence is coming only from the light degrees 
of freedom, which are at least tree-level improved 
through 0(a), one expects M = 0(a s AQCD), an 
order of magnitude smaller than what is found. 
This may be a renormalon shadow effect. 

6. CONCLUDING REMARK 

I have not made a unitarity triangle analysis 
with my world average lattice results. However, 
Fig. [ll] shows the result of an analysis |58j that 
employs lattice values similar to those quoted 



6 In the Formilab formalism, there are no true power law 
divergences (aside from the standard additive shift in the 
mass) as a — -> for fixed mg because the relativistic theory 
is regained. The coefficients (e.g., di) vanish fast enough 
to cancel the power divergence in the higher dim ensi onal 
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Figure 11. 

Allowed region for the vertex of the unitary tri- 
angle (CKM parameters p and fj) from Ref. |8§|]. 



here|] Clearly, our knowledge of the triangle is 
getting quite precise. Over the last dozen years, 
the allowed region for the vertex has shrunk in 
area by almost a factor of 20, as is shown dra- 
matically in Fig. 4 of Ref. |39| . We can be proud 
of the fact that lattice computations have played 
an important role in this accomplishment. 
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